Complex spatial organization in a simple model of resource allocation 



Damian H. Zanette 
Consejo Nacional de Investigaciones Cientificas y Tecnicas, 
Centra Atomico Bariloche and Insituto Balseiro, 8400 Bariloche, Rio Negro, Argentina 

(February 8, 2008) 



O 
O 
O 

(N 
•*-> 
O 

O 

CO 
CN 

o 

< 

d 



> 

m 
o 
o 

o 
o 



> 
X 



A dynamical model for the distribution of resources be- 
tween competing agents is studied. While global competi- 
tion leads to the accumulation of all the resources by a single 
agent, local competition allows for a wider resource distri- 
bution. Multiplicative processes give rise to almost-ordered 
spatial structures, thourgh the enhancement of random fluc- 
tuations. 

PACS numbers: 87.23. Ge, 87.23. Cc, 05.65.+b 



Though many real systems in the scopes of biology 
and the social sciences are well described by agent-based 
models with pair interactions Jl|^, qualitatively similar 
to the physical description of interacting particles, a large 
class of biological and social processes are driven by in- 
teractions mediated by an external actor, which generally 
bears no resemblance with the individual agents. This 
kind of processes — which, in connection with physics, 
may be assimilated to the evolution of globally cou- 
pled dynamical systems [Q — is typically found in systems 
where agents compete for resources. A particularly rel- 
evant instance of this situation is given by competing 
biological species. In this case, in fact, the interaction is 
rarely given by struggling events between individuals of 
different species, but rather by accessing simultaneously 
to the limited resources provided by the environment. 
The performance of each agent, i.e. of each species, is 
here measured by its ability to get resources, which de- 
pends both on the capabilities of each individual and 
on global features such as the total population of the 
species. The same scenario is found in some economical 
systems, for instance, in companies competing for finan- 
cial resources, usually administrated by banks, or even 
in scientific research projects competing for funds from 
a government agency. A key ingredient in the dynam- 
ics of these systems is that the ability of each agent to 
get resources at a given time depends, often strongly, on 
the resources assigned to the agent at previous stages. 
This can give rise to a multiplicative process that, in the 
absence of buffer mechanisms, leads to resource accumu- 
lation by a single agent. In the following, we consider a 
simple dynamical model that incorporates these elements 
and, in particular, study the effects of competition at lo- 
cal level. We find that local competition softens the pro- 
cess of resource accumulation and gives rise, through the 
enhancement of spatial fluctuations, to nontrivial struc- 
tures. 



Our system consists of an ensemble of N agents, each 
of them characterized by its productivity a,-. For fu- 
ture convenience, we consider the generic situation where 
all productivities are different. At each time step t, 
each agent is assigned an amount of resources rj(t), 
which is used to produce an amount of products Pi(t) = 
a.iTi(t). At the next time step, resources are distributed 
among agents in amounts proportional to their produc- 
tion, namely, 



n(t + i) = 



Pi(t) 



R(t+1) 



airi{t) 



R(t + 1), 



(1) 



where R(t + 1) are the total resources available at time 
t + 1. Model (|]) can be fully solved for arbitrary produc- 
tivities and initial conditions rj(0). In the first place, we 
note that rescaling r.i(t)/R(t) — » ri(t), Eq. (|l|) becomes 
independent of the total resources. We thus fix R(t) = 1 
for all t, so that J2j r j(t) = 1- hi this situation, the 
solution to Eq. (fy) reads 
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For asymptotically large times, ri(t) — > for all i, ex- 
cept for the agent with the maximal productivity, a max = 
maxijai}, which receives all the available resources. Due 
to the multiplicative effect of resource allocation accord- 
ing to production, all resources are in the long run ac- 
cumulated by the agent with the maximal productivity, 
giving rise to a sort of winner-takes-all state || . By anal- 
ogy with population dynamics we say that the remaining 
agents become extinct. In fact, in connection with bio- 
logical populations, this result is a realization of a well- 
known principle of ecology, namely, the principle of com- 
petitive exclusion |5]|l|] : in a system of biological species 
competing for the same resources, only one survives and 
the others undergo extinction. At moderately large t the 
resources assigned to each agent are well approximated 
by n(t) = (a i /a max ) t r i (0)/r max (0). Assuming that the 
productivities cti are drawn at random from a distribu- 
tion P a (a) and that, for simplicity, resources are evenly 
assigned at the beginning, rj(0) = A^ -1 for all i, the prob- 
ability distribution for the individual resources at a given 
time is 



P r (n) 



fit 



„l/t 



(3) 
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The dependence of this function on through the power 



is very weak for large i, so that on a wide interval 



of the variable we find P r (ri) oc i\ . For long times, 
thus, resources are distributed among agents following a 
power law with exponent — 1. Compare this result with 
Pareto's law of wealth distribution J(||]]. 

Model admits several variations, that may be of 
interest in connection with the description of some real 
systems. For instance, the extinction of all but one 
agent can be avoided by modifying slightly the alloca- 
tion method. If, at each time step, a small fraction p of 
the total resources is evenly distributed among the agents 
whereas the remaining is assigned according to produc- 
tion as above, the individual resources are always larger 
than p/N. In this situation, P r (ri) becomes stationary 
for long times and, if p is small enough, the power-law 
dependence quoted above is found again for Vj, > p/N . 

A second variation consists in assuming that the pro- 
ductivity can in turn depend on the individual resources, 
for instance, as ai(n) — a^A{r,i). If A{r) increases with 
r, the winner-takes-all effect is enhanced. The opposite 
case, where A{r) is a decreasing function, is more in- 
teresting. Lower productivities for higher resources may 
be the consequence of size effects, crowding, or loss of 
efficiency due to lower competition. In this case, the sys- 
tem evolves to a stationary situation where the agents 
whose productivity weights a? are above a certain thresh- 
old o^ in receive nonzero resources, whereas the agents 
with a° < a° nin become extinct. In the stationary situa- 
tion, all the surviving agents have the same productivity 
a. The third instance, where A(r) is a nonmonotonous 
function, can give origin to multiple stationary nontrivial 
solutions. 

In the following we focus the attention on a variation of 
model (|j) that introduces a spatial distribution of agents. 
Though the total resources are still allocated among all 
the ensemble, the agents compete at a local level only. 
The set of agents that compete with a given agent i de- 
fines its neighborhood Mi. Resources are assigned ac- 
cording to 



n(t + i) = 



1 



where 



m = £ 



in(t) 



(4) 



(5) 



is a normalization factor that insures that Y]^ rj (t) = 1 
for all t. We have performed numerical simulations of this 
system in ensembles of TV = 10 3 to 10 4 agents distributed 
on various geometries. The productivities a.; were drawn 
at random from a uniform distribution in (0,1). Note 
that, in Eqs. (|l|) and (^), a homogeneous rescaling of 
all productivities leaves the models invariant. Several 
distributions of initial conditions were tested, but there 



are no essential differences with the case where, at the 
first step, resources are evenly allocated, ri(0) = iV -1 . 
Therefore, we concentrate our simulations on this simple 
case. 

In the first place, we consider a one-dimensional ar- 
ray of agents with periodic boundary conditions. The 
neighborhood of the ith agent consists of its two nearest 
neighbors. We find that, for long times, approximately 
half of the agents become extinct and resources are evenly 
allocated among the surviving agents. Practically every- 
where along the array agents are ordered in an alternat- 
ing sequence of extinct agents and survivors. Occasion- 
ally, one finds two neighbor sites where both agents are 
extinct, but two survivors are never contiguous. In other 
words, in the neighborhood of a surviving agent no other 
survivor can be found [cf. model (@)]. The surprising 
feature of this asymptotic distribution is that the strong 
correlation between survival and productivity found for 
Eq. (|l|) is apparently lost in Eq. (|J). In fact, produc- 
tivities are distributed completely at random along the 
array, whereas the resulting distribution of survivors is 
highly ordered. At first glance, no connection can be es- 
tablished, for instance, between survivor sites and local 
maxima of the productivity or any spatial pattern caused 
by fluctuations in the distribution of cv^. 




FIG. 1. Asymptotic state on a 100 x 100-site lattice with 
nearest-neighbor competition and periodic boundary condi- 
tions. Full squares represent survivors. 

The same feature is found in a two-dimensional array 
with periodic boundary conditions, where the neighbor- 
hood of each site consists of the four nearest neighbors. 
Figure |l| displays the asymptotic state on a 100 x 100- 
site lattice (N — 10 4 ), where the full squares represent 
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survivor sites. Regular domains where survivors and 
extinct agents alternate in both dimensions, separated 
by worm-like boundaries formed by extinct agents, are 
apparent. We have also verified that different defini- 
tions of the neighborhood of a site, both in one and in 
two-dimensional lattices, produce similar almost-periodic 
asymptotic structures. Again, resources are evenly dis- 
tributed between the surviving agents and, in neighbor- 
hood of a survivor, no other survivor can be found. For 
instance, in a one-dimensional array where the neighbor- 
hood consists of four sites (nearest and next-to-nearest 
neighbors), the resulting structure is a periodic sequence 
of two extinct agents and one survivor. Occasional de- 
fects, with larger zones of extinction, are also found. 

In order to trace the origin of the almost-periodic struc- 
tures emerging at asymptotically long times we have 
inspected the successive states of the system from the 
earliest evolution stages, for the specific case of a one- 
dimensional array of agents with two neighbors per site 
and periodic boundary conditions. As above, productivi- 
ties are drawn at random from a uniform distribution and 
ri(0) = TV -1 for all i. We find that even at the first step 
of the evolution the distribution of resources already ex- 
hibits an almost-periodic sequence of relatively high and 
low values, in spite of the fact that productivities are spa- 
tially uncorrelated. Namely, for a substantial fraction of 
agents, we have [r i+ i(l)-r i (l)][r i (l)-r l _i(l)] < 0. This 
observation reveals that the simple dynamics of model 
is able to introduce strong spatial correlations to an 
initially uncorrelated state, even at the level of a single 
evolution step. 

To quantify this effect, we take a set of five uncorre- 
lated random numbers a±, ...,0:5 drawn from the pro- 
ductivity distribution, and consider the combinations 



Oi-i + on + a l+ i 

for i = 2,3,4. The quantities <Xj are then proportional 
to the resources received at t = 1 by three consecutive 
agents in the array [cf. Eq. ([|)]. Numerical realiza- 
tions of these quantities, over 10 6 independent choices of 
ax, ■ ■ ■ , c*5, show that (174 — <J^){a^ — 02) < with prob- 
ability pi = 0.762. Consequently, slightly more than 76 
% of the agents are expected to belong to a periodic se- 
quence of alternating high and low resources at the first 
time step. The remaining 24 % stands for defects in the 
periodic structure. The same kind of analysis can be per- 
formed for successive steps in the evolution. It is found 
that correlations are further enhanced by the dynamics. 
At the second evolution step, for instance, the probabil- 
ity for an agent to belong to the periodic structure grows 
to p 2 = 0.782. For later stages, we find p$ = 0.816, 
P10 = 0.839, p 20 = 0.854, and p 50 = 0.864. These proba- 
bilities saturate at p^ « 0.865. 
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FIG. 2. Survival probability as a function of the productiv- 
ity a for one dimensional arrays of 1000 agents and different 
neighborhoods. Dots correspond to average measurements 
over 100 realizations. The line shows the result of an inde- 
pendent calculation method (see text), for the case of two 
neighbors. 

While the emergence of almost-periodic structures 
seems to indicate that no correlations subsist in model 
([j]) between survival and productivity, the inspection of 
the evolution equations as well as of Eq. (^) suggests 
that some remaining correlation should however exist. 
In fact, low productivities still imply few resources, so 
that the minima in the periodic structure should prefer- 
ably coincide with sites with small aj. Moreover, de- 
fects in the periodic structures — where extinct agents 
are found in consecutive sites — should also correspond to 
low-productivity zones. Such remaining correlation can 
be characterized by the probability of survival as a func- 
tion of the productivity, p s (a). This probability is de- 
fined as the fraction of agents with productivity a which 
survive at asymptotically long times. We have numeri- 
cally measured p s (a), in series of 100 independent realiza- 
tions of the productivity distribution for ensembles with 
N = 1000, for various geometries and neighborhoods. 
Figure || shows the survival probability as a function of 
a for one-dimensional arrays with different numbers of 
neighbors. Here, the correlation between survival prob- 
ability and productivity is apparent. Agents with larger 
productivities are more likely to survive that those with 
small a. Note, however, that for a = 1 the survival prob- 
ability is less than unity, so that even with the maximal 
productivity there are chances of undergoing extinction. 
Conversely, for a — > (but a ^ 0), the survival prob- 
ability is finite. As the number of neighbors grows the 
probability is more concentrated towards larger produc- 
tivities, as expected. In the limit where the neighborhood 
extends to the whole array the original model (|l|) is re- 
covered, and the survival probability must vanish except 
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for a = 1. 

The survival probability can be calculated indepen- 
dently of the numerical realization of the model, using 
expressions as in Eq. (|^) . For a fixed value of and dif- 
ferent random choices of ot\, a%, and 05, an estimate 
at the first evolution step of the survival probability of 
the agent with productivity a% is given by the fraction 
of instances for which 173 > <J2, 04. Generalizing this pro- 
cedure for successive time steps, the estimate can be im- 
proved by considering later stages in the evolution. The 
line in Fig. |^ corresponds to this estimate at t = 50. It 
shows an excellent agreement with the numerical realiza- 
tions. 



1.0 



0.8 



o 



C3 
> 

U 

3 



0.4 



0.2 



■ one dimension 
A two dimensions 
* random graph 



• jfT 




0.6 



1.0 



FIG. 3. Survival probability as a function of the produc- 
tivity a for different geometries with four neighbors per site. 



We may summarize our main results as follows. First, 
the strong, deterministic correlation between survival 
and maximal productivity that characterizes the model 
with global competition (|l|) is replaced, in the model with 
local competition (Q), by a much weaker, probabilistic- 
like correlation. Agents with high productivity are more 
likely to survive but, even with the maximal productivity, 
extinction cannot be completely discarded. On the other 
hand, agents with very low productivity have a chance of 
survival. This conclusion should be relevant to the possi- 
ble applications of the present models, both to economy 
and to biology. Second, the loss of the above mentioned 
deterministic correlation is accompanied by the forma- 
tion of an almost-regular spatial structure, with a peri- 
odic alternation of survivors and extinct agents. This 
structure is explained by the emergence of strong spatial 
correlations out of the fully uncorrelated distribution of 
productivities, due to the very action of the evolution 
rules. Nontrivial correlations in the conditional probabil- 
ities for certain combinations of uncorrelated variables — 
a rather simple but counterintuitive phenomenon — have 
recently been pointed out for series of random numbers 
. The present model illustrates the occurrence of the 
same kind of phenomenon in a spatially extended dynam- 
ical system. 

This work has been partially carried out at the Abdus 
Salam International Centre for Theoretical Physics (Tri- 
este, Italy). The author thanks the Centre for hospitality, 
and G. Abramson for useful comments. 



To compare now several geometries, we have measured 
the survival probability in systems where each agent 
has the same number of neighbors — specifically, four — 
on different kinds of arrays, namely, a one-dimensional 
array (nearest and next-to- nearest neighbors), a two- 
dimensional lattice (nearest neighbors) and a random 
graph with four connections per site. The results, shown 
in fig. [?] corresponds to averages over 100 realizations in 
systems with N = 1000 for one dimension and random 
graphs. In this latter case, the graph topology is chosen 
anew at each realization. For two dimensions, the re- 
sults correspond to 10 realizations on a 100 x 100-array. 
Though the general trend of p s (a) is qualitatively the 
same for the three cases, some systematic differences are 
apparent. For instance, the survival probability in one 
dimension is appreciably lower than in the other geome- 
tries for low productivities. In this range, remarkably, 
the data for random graphs lies between those for one 
and two dimensions. For a > 0.8, instead, the values 
of p s (a) are hardly distinguishable within our numerical 
precision. 
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